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ABSTRACT 

We examine the use of the CMB's TE cross correlation power spectrum as a com- 
plementary test to detect primordial gravitational waves (PGWs). The first method 
used is based on the determination of the lowest multipole, £o, where the TE power 
spectrum, Cf^, first changes sign. The second method uses Wiener filtering on the 
CMB TE data to remove the density perturbations contribution to the TE power 
spectrum. In principle this leaves only the contribution of PGWs. We examine two 
toy experiments (one ideal and another more realistic) to see their ability to constrain 
PGWs using the TE power spectrum alone. We found that an ideal experiment, one 
limited only by cosmic variance, can detect PGWs with a ratio of tensor to scalar 
metric perturbation power spectra r = 0.3 at 99.9% confidence level using only the 
TE correlation. This value is comparable with current constraints obtained by WMAP 
based on the 2a upper limits to the B-mode amplitude. We demonstrate that to mea- 
sure PGWs by their contribution to the TE cross correlation power spectrum in a 
realistic ground based experiment when real instrumental noise is taken into account, 
the tensor-to-scalar ratio, r, should be approximately three times larger. 

Key words: cosmic microwave background - polarization - gravitational waves - 
cosmological parameters 



1 INTRODUCTION 

Primordial gravitational waves (PGWs) polarize t he cosmic 
' microwav e back g round (CMB) (s ee for example Basko & 
Polnarcvl (|l980|); |PoInarcv| (tlQSSl); ICrittenden et al.| ([19931) 



Frcwin 




IColcs ct al, (1990; Kamionkowski 



dent components: E-mode (grad) polarization and B-mode 
(curl) polarization. B-mode pol arization can o nly be gener- 
ated by PGWs ( see for example ISeljakl ([L997|) ; feeliak fc ZaL 
darriaga (1997); Kamionkowski & Kosowsky (1998)), there 



ISeljak fc Zaldarriagal (1997); 



Kamionkowski & Kosowsky ( |1998 )7 ^askaran et al. (2006); 



fore most CMB polarization experiments which are search- 
ing for evidence of PGWs focus on measuring the BB power 
spectrum. However the TE cross correlatio n power spectrum 



Keating ct al.| (2006)). Current experiments are using the 
polarizat ion of the CMB to search for this PGW back- 
ground (iTaylor et al.| ( |2004| ) 



offers another method to detect PGWs (Crittenden et al 
( 1995| )). The TE power spectrum is two orders of magni 



Bowden et al 



(2004); yoon 



et al. (2006). This polarization can be used as a direct test 
of inflation. An alternative probe of the inflationary epoch 
which does not use the PGW background was studied by 
(Spergel & Zaidarriaga (1997)). Thi s probe was u sed in re- 



tude larger than the BB power spectrum and it was sug- 
gested that it may therefore be easie r to detect gravi t ationa l 



waves in t h e TE power spectrum (Baskaran et al 
Grishchukl (^007|)). 



(2006); 



cent analyses by the WMAP team (Peiris et al. (2003)) to 



provide plausibility for the inflationary paradig m. This pa- 
per presents a test similar in spirit to that of Spergel & 



Zaidarriaga) ( |l997| ). 

CMB polarization can be separated into two indepen- 
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In this paper we first discuss the method of detection 
of PGWs by measuring the TE po wer spectrum for low £ . 
This method, originally proposed in Baskaran et al. (2006), 
is based on a measurement of lo, the multipole where the 
TE power spectrum first changes sign. Hereafter we will call 
this method "the zero multipole method". The TE power 
spectrum due to density perturbations is positive on large 
scales, corresponding to ^ < ^o, changes sign at £ = £o, and 
then oscillates for £ > £o, while for PGWs the TE power 
spectrum must be negative for small £ and then also os- 
cillates for larger £. The current best set of cosmological 
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parameters, obtained in Spergel et al. (2007), gives, in the 
absence of PGWs, lo ~ 53. Therefore, the measurement of 
the difference between the multipole number, £o, where the 
TE power spectrum changes sign, and ^ = 53 is the way to 
detect PGWs. We will then consider an alternative method 
based on Wiener filtering, removing the contribution to the 
TE power spectrum due to density perturbations. Since the 
TE power spectrum due to PGWs is megative on large scales 
a test of negativity of the resulting TE power spectrum is a 
test of PGWs. In this paper, we present an analysis of both 
of these methods, based on Monte Carlo simulations. 



At the present time, the main priority and the main 
challenge in CMB polarization observations is the detection 
of the PGW background via the BB power spectrum. In con- 
nection with BB experiments, the methods based on the TE 
cross correlation can be considered as very useful auxiliary 
measurements of PGWs because systematic effects in TE 
measurements are not degenerate with those in BB measure- 
ments. For example, T/B leakage or even E/B leakage could 
swamp a detection of BB, w hereas T/E leakage w ould be 
small and well controlled (see ^himon et al (2007)). These 
BB systematics could falsely indicate a detection of PGWs, 
but measurements of the TE power spectrum provide insur- 
ance against such a spurious detection. Additionally, galac- 
tic foreground contamination affects BB and TE in different 
ways, which enables us to perform powerful cross-checks and 
subtraction of foregrounds in BB measurements. 



Another advantage of TE measurements for experi- 
ments which measure a small fraction of the sky, is related 
to the fact that a significant contaminant to the B modes 
is caused by E/B mixing. This limits the power s p ectru m 
of PGWs that can be detected ( |Challinor fc Chon| ( |2005| )). 
The E- modes are practicly unaffected by E/B mixing so, 
in contrast to the BB measurements, the TE power spec- 
trum should be nearly the same for both full and partial sky 
measurements. 

The plan of this paper is the following. In Section we 
introduce the primordial power s pect ra of scalar (de nsity) 



and te nsor (PGW) pertu rbations ( 2.1 ). Then fo llowing Crit- 
tenden et al. (1995) and Baskaran et al. (2006), we explain 
why the sign of the TE power spectra for scalar and tensor 
perturbations is opposite for large scales (2.2). In Section ^ 
we describe in more detail the zero multipole method for the 
detection of PGWs. In Section ^, we describe the method for 
detection of PGWs based on Wiener filtering along with the 
statistical tests used and a comparison of the tests. In Sec- 
tion ^ we present results of numerical Monte Carlo simula- 
tions for two toy experiments. In the first toy experiment we 
neglect instrumental noi se a nd the uncertainties are limited 
only by cosmic variance (5.1). In the second toy experiment. 



along with cosmic variance, we take into account instrumen- 
tal noise whi ch is comparable to real noise in current ground 
experiments (5.2). For comparison, we also present resul ts o f 



simulations for the two satellite experiments, WMAP (5.3) 
and Planck (5.4). In Section ^, we compare the the signal- 
to-noise ratio of the TE measurements with those of BB 
measurements. 



2 TE CROSS CORRELATION 

The power spectrum of TE correlations is determined by 
primordial power spectra of scalar and tensor perturbations 
and time evolution of these perturbations during the epoch 
of recombination. 



2.1 Primordial Power Spectra 

The primordial power spectra describing the initial scalar 
(density) perturbations (denoted by s) and tensor (PGW) 



perturbations (denoted by t) are (see, for example, Spergel 
(|2007|)) 



et al. 



Ps{k) = h- 



Us — l+^OLs log(fc/fco) 



(1) 



where fco = 0.002 Mpc ^, this value of ko is obtained by 



fitting of CMB data ( [Smith et al.| ( |2006| )). The variables Us 
and rit are the scalar and tensor spectral indices respectively. 
The variable as is the running of the scalar spectral index. 
In terms of As and At, the tensor-to-scalar ratio, r, is 



As 



Ptjko) 
Ps{ko) 



(2) 



The location of £o is determined by the parameters nt 
and r. In this paper, we do not specify particular cosmologi- 
cal models considering the generation of primordial spectra, 
Ps (k) and Pt (fc) , which means that for our purposes we con- 
sider Us, nt, and r as independent parameters. (This is not 
true if we use some particular cosmological model. For ex- 
ample, in standard inflation models, the parameters nt and 
r are rela ted by the consisten cv relation, nt — — r/8 (see, for 



example, Peiris et al, (2003)).) In other words, we consider 



all parameters n^, nt, and r as independent except in 5.1 
and where along with model-independent we give also 
model-dependent constraints on r. 



2.2 Opposite Signs of Scalar and Tensor 
Perturbations to TE Correlation 

Taking into account that scalar and tensor perturbations are 
not correlated, the TE power spectrum is simply a sum of 
two TE power spectra for scalar and tensor perturbations 
correspondingly. 

First, the physical motivation for the difference in the 
cross correlation contributions produced by scalar and ten- 
sor perturbations for small £ was demonstrated and physi- 
cally interpreted for the cross correla t ion o f the Stokes pa- 
rameters T and Q in Crittenden et al. (1995). For scalar per- 
turbations the Stokes parameter Q contains only E-modes, 
hence the TE correlation is identical with the TQ correla- 
tio n and is positi v e for small £. As was then emphasized 
in Baskaran et al. (2006), the sign of the TE correlation for 
tensor perturbations is negative for small £. The simple qual- 
itative physical interpretation of the fact that the contribu- 
tions of the TE correlation are different for scalar and tensor 
perturbations is the following. For both scalar and tensor 
perturbations, the temperature fluctuations, T{£), for small 
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I (when oscillations of T{t) are absent) are proportional to 
the metric perturbations h at the moment of recombination, 
while the E-mode fluctuations, E{t), are proportional to h 
at the moment of recombination. Hence, the TE correla- 



tion is proportional to hh oc 



dri 



Taking into account the 



growth of scalar perturbations and tensor perturbations de- 
cay, one can see that the contributions to the TE correlation 
for scalar and tensor perturbations are opposite 



To und erstand this in more detail, following Baskaran 



(^006), we consider the multipole expansion of the TE 



et al, ^ 

cross correlation with coefficients Cj^ . These coefficients 
are related to the spherical harmonic expansion coefficients 
of the temperature anisotropy and polarization by 



the TE cross correlation power spectrum changes sign as a 
function of I. 

Thus the TE cross correlation, due to density pertrba- 
tions, must be positive at lower I (as mentioned above, the 
TE cross correlation in absence of PGWs changes sign at 
to ^ 53). If we were able to separate them we could use 
this signature for detection of PGWs. However, even with- 
out such separation the presence of PGWs manifests itself 
in the value of lo, which is the smallest i where the total TE 
correlation power spectrum (scalar plus tensor) changes its 
sign. Thus, the sign of the TE correlation is a very promi- 
nent signature of PGWs. For this reason, in the next section, 
we investigate the dependance of on r, Us, and nj. 



C 



7 



(3) 



where the brackets denote averaging over all possible statis- 
tical realizations. The statistical properties of the CMB field 
in general, and the TE cross correlation specifically, follow 
from the statistical properties of the underlying scalar or 
tensor metric field. Assuming gaussianity together with sta- 
tistical isotropy and homogeneity, the TE cross correlation 
takes the form 



C 



aT,e{k)aE,i{k) 



(4) 



where aT,i{k) is the contribution from temperature pertur- 
bation while aE,t{k) is the contribution from E-polarization. 
The integration over k takes into account the contribution 
from all the possible wavenumbers. 
It was shown in Baskaran et al. 



(2006) 



aT,i{n) ~ hk[ri) 



a-E, 



tin) 



dhk{rj) 



drj 



(5) 



(6) 



where {rj) is the mode function of the metric perturbation, 
and rjrcc is the conformal time at recombination. It follows 
that the TE correlation is approximately 



Ci oc 



dnF{e., k) 



dr] 



(7) 



where F{£, k) is a strictly positive function which peaks at 
~ f « k {ritoday —rjrec)- Hcuristically, the function F{£,k) 
projects the fc-space onto the ^-space. Therefore the sign 
of the integral in the RHS of Eq. ^ evaluated at around 
£ ~ k {ritoday — Tjrec) determines the sign of Cj^ on large 
scales. 

The adiabatic decrease of the gravitational wave am- 
plitude upon entering the Hubble radius is preceded by the 
monotonic descrease of the gravitational wave mode func- 
tion hkijf) as a function of 77. Since hkijf) is decreasing the 
integral on the RHS of Eq. 0is negative. The RHS of Eq. 0is 
negative for k {rjtoday — rjrec) < 90 since hk is decreasing over 
that range. Therefore, for £ < 90 the correlation Cf^ must 
be negative. For larger ^s, the F{£, k) in Eq. ^ and, hence. 



3 DEPENDANCE OF £0 ON PARAMETERS OF 
PGW POWER SPECTRUM 

The method of detecting PGWs which implies a calculation 
of £q, where the TE power spectrum first goes to zero, will 
be called hereafter as the zero multipole method. We take 
into account uncertainties in determination of Ces which are 
unavoidable in any experiment: 



i2£+l)fsky 



{crf + 



icr+NDicr+ND 



(8) 



(see, for example, Dodelson (2003)). Even in an ideal ex- 
periment, when we neglect instrumental noise {Ni = 0) and 
measure the full sky {fsky = 1), we still have uncertain- 
ties related with cosmic variance (which arises from the fact 
that we have only one realiza tion of th e sky in CMB mea- 
surements) (see, for example, Dodelsor (2003)). For a more 



realistic experiment, we take into account noise and partial 
sky coverage (see Section ^). Over small multipole bands it 
is reasonable to approximate the power spectrum as linear. 
In the range 20 ^ £ ^ 70, it seems reasonable to use a linear 
approximation for {£ + 1)Ci/2tt. It seems unlikely that in 
this range any deviations from a linear approximation can 
be larger than mentioned above uncertianties. 

Plots of Cj^ for different values of r are shown in Fig. 
|l| plotted for nt = 0. It can be seen that a linear fit to the 
TE power spectrum do well approximates {£ + l)Cj^ /2-k 
near £(,. 

Thus near £q, {£ + l)Cf^ /2it can be approximated as 
a line with negative slope a — b£, where a and b are positive 
real numbers. For any set of experimental data, we can find 
a and b by applying a least squares fit. The values a and 
6 corresponding to the best fit obviously can be used for 
prediction of £0 = a/b. This value, £0, can then be used to 
constrain the parameter r under some assumptions about 
spectral indices Us and rit- 

We need to investigate how £0 depends on the cosmolog- 
ical parameters r, nt, Us, As — Ps{ko), and the optical depth 
to reionization, r. T he value of fp for a st andard ACDM cos- 
mology described in 3p_ergel et al. (2007) as a function of nt 
and r is shown in Fig. ^ and ^ AH power spectra were gen- 
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Figure 2. Plot of crossover multipole number £o for different 
values of nj. nt = —0.5 to 0.5 with spacings of 0.1. The dashed 
red line correspond to nt = 0. 



erated with the code CAMbQ ( [Lewis et al] ( |2000| )). If r = 0, 
£o = 53, while if r = 0.3 (the WMAP3 upper limit on the 
tensor-to-scalar ratio) and nt — 0, we find that £o = 49. 

From Fig. ^ and ^ one can see that £o decreases with 
increase of r. This effect is more pronounced for smaller nt- 
For example, if r = 0.3, then £o = 52 for nt = -1-0.5, £o = 49 
for nt = 0, and £o = 38 for nt — —0.5. The fact that £ is 
discrete (the plots are composed of a set of step functions) 
puts limitations on using this method for determination of r. 
For example, a value of r = 0.01 and r = will most likely 
correspond to the same £o and therefore no matter what 
the sensitivity is this method cannot distinguish between 
absence of PGWs and PGWs corresponding to such small 
r. For nt = (the Harrison-Zel'dovich scale-free spectrum), 
S£o = — 1 corresponds to Sr of 0.08. For negative values of 



see http://camb.info on web 



nt, S£o = — 1 corresponds to smaller Sr. For nt — —0.5, for 
example, 5£o = — 1 requires Sr ~ 0.02. 

The effect of variations of the scalar spectral index ris 
on £o is opposite: A decrease of £o with increase of r is more 
pronounced for larger ns, however we do not need to worry 
about it because ris is well constrained by the observations 



of TT and EE power spectra (see for example jpergel et al 
(^200?!) ) along with Ly-a measurements (see for example 
et al.| (2004)). Thus everywhere in this paper we us e n. 



Vie] 



0.95 (the value given by WMAP3 [Spergel et al.| ( [2007] )) with 
no running of the scalar spectral index, a^. A change of 0.2 
in the running of the scalar spectral index has no effect in 
the value of £o when using fco = 0.002 Mpc"^ Even if it did 
have an effect on the value of £o, it can be well constrained 
by the same observations that constrain ns . 

The value of £o does not depend on As, because for 
fixed r. At must change by the same factor as As leaving £o 
unchanged, i.e. any rescaling of the primordial power spec- 
tra does not change £q. The same thing happens when one 
varies the optical depth to reionization, r. The value of £o 
is in the range where the TE power spectrum for scalar and 
tensor perturbations depend on r in the same way for in- 
stantaneous reionization (they are damped by the factor, 
exp(— 2r), since the relevant scales are wi thin the cosmolo g- 
ical horizon at the time of reionization (Dodelson (2003)). 
Thus any variation of r can be considered just as a rescaling 
of the TE power spectrum which leaves the value of £o un- 
changed. It is possible for differen t reionization hist o ries t o 
cause a change in £o as shown in Kaplinghat et al, (2003), 
however we will assume instantaneous reionization for the 
purpose of this paper. 

Thus, even if we cannot separate the contributions of 
scalar and tensor perturbations to the TE power spectra, 
PGWs still leave their imprint on the value of £o. In the next 
section, we will consider the possibility of such separation 
with the help of Wiener filtering. 
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4 WIENER FILTERING OF THE TE CROSS 
CORRELATION POWER SPECTRUM 

Wiener filtering has been used often in the case of CMB 
data analysis. For example, it was used to combine multi- 
frequency data in order to remove foregro unds and extract 



the CMB signal from the observed data (Tegmark & Efs- 



tathiou (1996); Bouchet et al. (1999)). Here we examine the 



use of the Wiener filter to subtract the PGW signal from the 
total TE correlation signal. This is done because the Wiener 
filter reduces the contribution of noise in a total signal by 
com parison with an estimation of the desired noiseless sig- 
nal (Vaseghi ( 200(j )). In our case, the signal is the one due 
to PGWs only, and the signal contributed by density per- 
turbations is considered to be "noise" . 

The observed signal can be written as 



C 



{a*E 



) 



(9) 



where s and t refer to the contributions to the power spec- 
trum due to scalar and tensor perturbations respectively. 
The values aE,em and ax^em refer to the spherical harmonic 
coefficients of the temperature and polarization maps. In our 
application to TE correlation, we consider the Wiener filter, 

WTE,t'- 



Wte,! 



C 



\C, 



i.t 



(10) 



The filtered signal, ax,«m {^or X — T and E), is obtained 
from the measured signal, ax, em, as 



o-x,em 



,1/2 
' TE.i 



(11) 



In this paper, we assume the Wiener filter is perfect, in the 
sense that it leaves the signal due to PGWs only. We then 

-yTE 



get, for the filtered multipoles Cf ^ 



ryTE _ / '* ' \ 

'^ejat — \0'T,lmO.E,tm) 



= Wte ^iC^i^ = 



(~,TE 



CJ 



I — ^i. 



(12) 



In practice this is not true, because we are trying to 
determine Cff, which is not known in advance. Neverthe- 
less, the assumption that the Wiener filter is perfect is good 
as a first approximation and illustrates the detectability of 
PGWs with the help of TE correlation measurements. 

The filtering can reduce the measured signal to the de- 
sired signal, but, since we are trying to remove the density 
perturbations and not the actual noise, we can not reduce 
the measurement uncertainties. These uncertainties in Cf'^ 
are then entirely determined by the noise in the original 
signal. 

We have shown that the TE power spectrum due to 
PGWs is negative on large scales, hence a test determining 
whether the Wiener filtered power spectrum is negative or 
not is a probe of PGWs. 

There are three different statistical tests we use to see 
if we can measure a negative TE power spectrum. The first 
test is a Monte Carlo simulation to determine signal-to-noise 
ratio, 5'/A'^ (Section [l.l[). The other two tests are standard 



non-parametric statistical tests; the sign t e st (Se ction 4.2) 
and the Wilcoxon rank sum test (Wilcoxon ( 194£| )) (Section 

For all of our tests, we calculate a random variable. If 
the data satisfies the hypothesis that r = 0, we can calcu- 
late the mean and uncertainty in the variables. If we make 
one realization of data, the random variable is determined 
from its distribution. Because we are not using any real ob- 
servational data, we must run a Monte Carlo simulation to 
reduce the risk of randomly getting a value for the variable 
taken from the outlying area of its distribution. To do this. 



the filtered multipoles, C, 



TE 

e,fiit. 



are randomly chosen from a 



gaussian distribution with mean Cf f and standard devia- 
tion ACT'^ where 



(ACj 



(2£+l)fsky 

{Cr + NDiCr + Nf) I (13) 



(see, for example, Dodelson (2003)), the variable faky refers 
to the fraction of the sky covered by observations and A''^ is 
the effect i ve pow er spectrum of the instrumental noise (see 



Dodelson 



(2003) for details on how A'^ is related to actual 



instrumental noise). 

Our determination of Cjf is dependent on £. However, 
for two of our tests we ignore the value of £ in the calcula- 
tion of the random variable. We assume that the calculated 
random variable is gaussian. In order for this to work, the 
random variable must be calculated from gaussian variables. 
The errors on the multipoles for the "ideal" toy experiment 
are large enough so that we can assume the multipoles are 
taken from a single distribution and not from a distribution 
that depends on £. 



4.1 Monte Carlo S/N Test 

For this test, the random variable we calculate, S/N, is de- 
fined as 



C 



(14) 



The reason why the sum in this equation is taken in 
the range 2 < / < 53 is because only in this range 
S(7n(Cj^ (scalar)) = —sgn(Cj^ (tensor)). In other words, 
if we include higher multipoles we confront with a danger of 
a false detection, because the total TE power spectrum is 
negative for ^ > 53. 

The value of S/N is gaussian distributed because it is a 
sum of many modes of squares of gaussian distributed val- 
ues, Ce = a| We approximate each Cf^ as being gaussian 
distributed for the purpose of this paper. For each set of pa- 
rameters we run this simulation one million times to deter- 
mine the mean, {S/N}, and standard deviation, (Ts/n- The 
mean of this distribution is determined by the preassumed 
value of r, while the standard deviation is determined by 
parameters of the experiment and gives the confidence level 
of detection. We run such Monte Carlo simulations for dif- 
ferent values of r to determine in what range of r we can 
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detect PGWs. When then using real observational data, we 
can compare the actual value of S/N with the results of 
Monte Carlo simulations to infer the likelihood, as function 
of r, which determines the probability that r 7^ 0, or that 
PGWs exist at detectable levels. 



4.2 Sign Test 

The sign test is a test of compatability of observational data 
with the hypothesis that r = 0. If we do have r = 0, then 
Cf fin will be equally distributed around zero. Application 
of this test to the filtered data is very simple. In practice, 
all observational data are distributed between several bins 
and the averaging of the signal is produced in each bin sep- 
arately. Let Nbins be the number of such bins. The sign test 
actually gives the probability that in N- bins the average 
is negative and in N+ = Ntins — N- it is positive, if r = 0. 
This probability, P, is given by the binomial distribution 



^^^^lUlLo.S^'— (15) 



The probability that the hypothesis r = is wrong is 



P(r-/0)R.l-2^P(i 



(16) 



The value X^jJq P{i) is the probability that we would get 
S5 A'^+ positive values given r — 0. This is the same as the 
probability of getting ^ negative values given r = 0. 
Therefore our confidence that r 7^ is just 100% minus 
the sum of the probabilities describe above (the probability 
that the A''-|- is closer to the mean, Ni,ins/'2, if r = 0). This 
equation only makes sense if A'^+ < Nbins /i, since that is 
required for r > 0. If A'^+ > Ntins/'2, that would imply 
r < 0, which is not physical. We would have to interpret the 
result as a random realization of r J5 0, with the most likely 
result of r = 0. Therefore we would not be able to say r 7^ 
with any confidence. 

Let us consider the following example: we put all mea- 
surements of Cf^ into 11 bins and in three of them the 
average is positive. In this example, the probability that the 
hypothesis r = is wrong is equal to 89%. 

One possible drawback of this method is that it does not 
take into account any measure of the signal-to-noise ratio of 
individual measurements. As we show in Section i.4, it is 



possible to have two completely different sets of data with 
the same probability of having r = 0. This test is also unable 
to make any prediction as to the value of r, only that it 
differs from zero. 



4.3 Wilcoxon Rank Sum Test 

This statistical test deals with two sets of data. The first set 
of data is taken from a real experiment which measures Cf^ 
with some unknown r. The second set of data is generated 
by Monte Carlo simulations (see Section O) with r = 0. 
The objective of the Wilcoxon rank sum test is t o give the 
proba bility that the hypothesis r = is wrong (Wilcoxon 
( [l945|) ). 



First, we choose some random variable U, whose proba- 
bility distribution is known if r = 0. For that, let us combine 
all data from first set with ni multipoles and second set with 
722 multipoles into one large data set, which obviously con- 
tains ni +n2 multipoles. Then, we rank all multipoles in the 
large data set from 1 to ni +n2 according to their amplitude 
(rank 1 for the smallest and rank ni + 712 for the largest). 
Now, the variables Ri and R2 are defined as the sum of the 
ranks for the first original data set and the second original 
data set, correspondingly. Finally, the variable U, is 



U = min(C/i, f/2), where 
U^ = i?i - n,(n, + l)/2. 



i = l,2 (17) 



If all multipoles of the first data set are larger than all mul- 
tipoles of the second data set, then Ui = nin2 and U2 = 0. 
It is not difficult to show that Ui + U2 ~ nin2 . If both sets 
of measurements have no evidence for PGWs, (Ui) = {U2}- 
It is also simple to see that Ui + U2 = nin2- 

It is important to emphasize that the ranks of multi- 
poles are random variables because all multipoles themselves 
are random variables, hence Ui, U2, and U axe random vari- 
ables. If ni -I- 712 is large, the distribution of U can be ap- 
proximated as a gaussian with a known mean and standard 
deviation. In this approximation we have 



mu 



cru = 



7ll7l2/2 



711712 (ni +712 + 1) 

12 



(18) 
(19) 



In some cases, instead of U, the variables Ri or R2 are used. 
The reason U is used here is because mu is symmetric in the 
data sets. If r = in both sets of data, then the distributions 
of Ui and U2 are the same, no matter what 711 and 712 are. 
The distributions of Ri and R2 would be the same only if 
Til = 712 . The probability that the first data set corresponds 
to r 7^ obtained from the test in which Ri or R2 is used is 
the same as if U is used. 

Since this test requires Monte Carlo simulations for the 
second set of data, we ran this test many times for many 
different data sets to get an accurate mean value for U. 

To reject the hypothesis r = means to detect PGWs. 
Using the Wilcoxon rank sum test the allowable value of 
r is determined, if instead of comparing with simulated 
data with r = 0, we compare with simulated data with 
r = ro 7^ 0. In order to get a range of allowable values 
for r, we need to run multiple Monte Carlo simulations with 
multiple values for ro. This is where the assumption that the 
Cf f are from a random distribution that is independent of 
£ is used. This implies that the ranks are random variables. 
If the errors on the Cf f are small enough, then the ranks 
will be predetermined. Therefore, our assumption about the 
distribution of U will not be true and the test would have 
to be modified. Fortunately, this is not the case for even an 
experiment only limited by cosmic variance. 

To illustrate how this test works, let us consider the 
following example. Assume there are 4 multipoles in the first 
set of data and consider that r — 0.3 is the correct value. 
There are also 4 multipoles in the second set of data (which 
for sure corresponds to r = 0). All quantities below are 
expressed in fiK^ . The value for the first data set are Ciif = 
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-0.005, Cii^ = 0.02, Cjo^ = -0.015, and Cjo^ = -0.01. 
The values for the second data set are Cio^ = 0.03, C^o^ = 
0.003, Cio^ = -0.02, and Cl^ = -0.003. A ranking of 
multipoles gives the ordering from lowest to highest, with 1 
referring to the first data set and 2 referring to the second 
data set, as 21112212. This results in i?i = 2 + 3 + 4 + 7 = 16, 
;7i = 16 - 10 = 6, and f/2 = 16 - 6 = 10. Therefore U = 
min{10,6} = 6. For ni = n2 = 4, to reject the hypothesis 
that r = at 95% conf idence lev e l, U\ should be less than 
one (see, for example, Lehmann ( 1975 )). In this example, 
since Ui = Q > 1, the first set of data cannot be considered 
as a detection of PGWs. 



4.4 Comparison of Tests 

The S/N test is greatly affected by outlying measurements. 
A measurement of one large negative multipole could falsely 
imply a detection. Both the sign test and the Wilcoxon rank 
sum test are not affected by individual outlying measure- 
ments. In the sign test, the value of individual measurements 
is irrelevant, because the test is sensitive only to the sign of 
individual measurements. The Wilcoxon rank sum test is af- 
fected by outliers, but considerably less than the S/N test. If 
the outlier is larger (or smaller) than every other multipole, 
its rank does not depend on its particular value. 

If we have two completely different sets of data, the 
main disadvantage of the sign test, as mentioned in Section 



4.2, is that it could give the same result, while for the two 
other tests the chance to obtain the same value of r is neg- 
ligible. For example, one set of data, consisting of 4 small 
negative multipoles and 4 large positive multipoles, gives the 
same result as another set of data, consisting of 4 large neg- 
ative multipoles and 4 small positive multipoles. The S/N 
test gives two very different values of S/N for these two sets 
of data. We can also use the Wilcoxon rank sum test to 
compare these two sets of data. In this case U — IQ = ^mu, 
which corresponds to a confidence level of hypothesis that 
r = of less than 10%. 

With observational data, the sign test can be applied 
and does not require any Monte Carlo simulations (which 
could be considered as an advantage of this test). The S/N 
test requires Monte Carlo simulations, but only for the dis- 
tribution of the random variable S/N. The Wilcoxon rank 
sum test requires large Monte Carlo simulations and com- 
bines the data sets generated by these simulations with ob- 
servational data. In other words, Monte Carlo simulations 
are absolutely necessary after obtaining observational data, 
which may be considered a disadvantage of this test. Thus, 
each of the three tests has advantages and disadvantages, 
suggesting that the best way to work out observational data 
is to apply all these three tests. 



5 DISCUSSION AND RESULTS 



Baskaran et al. (2006) used equal amplitudes of scalar and 



tensor perturbations to sharpen the discussion in their plots. 
They defined the tensor-to-scalar ratio, R, as the ratio of the 
temperature quadrupoles, R = C2J /C2J and set R = 1 



in this paper is approximately twice as large as their def- 
inition of R. The exact relationship between r and R will 
depend on the cosmological parameters used. This means 
that _R = 1 is equivalent to r ~ 2, which has currently been 
strongly ruled ou t by WMAP i n com bination with previ- 
ous experiments (Spergel et al. (2007)). We need to see if 



Using stand ard WMAP3 cosmological parameters ( Spergel 



et al. ( 2007 )), the definition of tensor-to-scalar ratio, r, used 



this method can detect a value of r that is currently within 
the limits. We assume that there is no foreground contam- 
ination. In reality foregrounds affect the measured location 
of £0 (we will consider the effects of foregrounds on £0 else- 
where) . For the experiments that do not observe the full sky, 
correlations between multipoles must be taken into account. 
The multipoles are binned together of such width that the 
correlations between the bins are sufficiently small. 

Two different toy experiments, along with the two satel- 
lite experiments WMAP and Planck, are considered to con- 
strain r. The first toy experiment is a full sky experiement. 
It is idealized in two aspects. The first idealization is that 
we can take measurements over the full sky while the second 
idealization is that we assume there is no detector noise. The 
only uncertainty is due to cosmic variance. Such experiment 
represents the best limit to which the gravitational waves 
can be detected by the CMB TE correlation. This toy ex- 
periment is close to a space-based experiment with access 
to the full sky. It is similar to what the Beyond Einstein in- 
fiation probe would be able to detect. This toy experiment 
will be hereafter referred to as the ideal experiment. The 
second toy experiment is a more realistic one. In this exper- 
iment, measurements are on 3% of the sky, the frequency 
is 100 GHz, and the duration of the experiment is 3 years. 
The noise in each detector of the 50 polarization sensitive 
bolometer pairs can be described by their noise equivalent 
temperature (NET) of 450 fiKy^. The detectors' beam pro- 
file is assumed to be gaussian and and it is described by the 
width at half of the maximum sensitivity, abbreviated as 
FWHMof0.85°. 

This second toy experiment is similar to current ground- 
based experiments and the constraints from this experiment 
represent those that can and will be obtained in the next 
several years. This will be referred to as the realistic exper- 
iment. 

The predicted errors for Planck are based on using the 
100 GHz, 143 GHz, and 217 GHz channel in the High Fre- 
quency Instrument (HFI) . The numbers are gotten from the 
Planck science case, the "bluebook"^. The WMAP noise was 
obtained by 3 years of observations of the Q-band, V-band, 
and W-band detectors. 



5.1 Ideal Experiment 

A plot of TE for r = 0.3 and nt — 0.0 with error bars for 
£ binned in bins of width A£ — 10 is shown in Fig. ^. This 
figure separately shows the contribution to the TE mode of 
density perturbations, contribution of PGWs with r — 0.3, 
when the TE power spectrum due to density perturbations 
is approximately five times larger than the power spectrum 
due to PGWsat£<£o. 

The Monte Carlo simulation for the calculation of £q, 
with an input model of r = 0.3 and nt = 0, results in the 



http: / /www.rssd.esa.int/index.php?project=Planck 
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Figure 4. The black line is the total TE mode with a r = 0.3. 
The red line is the contribution from PGWs only while the light 
blue line is the contribution from the density perturbations. Blue 
is the error bars for the ideal experiment binned in intervals of 
A£ = 10. 




0.0 0.1 0.2 0.3 0.4 0.5 
r 

Figure 6. The signal-to-noise ratio for the zero multipole method 
are shown as the solid black, for ideal experiment, and dashed 
blue, for realistic experiment, lines. The signal-to-noise ratio for 
realistic measurements of the BB power spectrum is shown as the 
dash-dot red curve. For all curves nt = 0. 




Figure 5. This is a plot of the allowed r and nt for the Icr region 
of £o for the ideal experiment. The white is the la region while 
the black is the forbidden region 



value of lo ~ 49 with an uncertainty of A^o ~ 1-3. A con- 
tour plot of the limits on the resulting measurement of r is 
shown in Fig. ^ The white is the allowed region for r and nt 
that falls within the Icr errors of lo. The black is the region 
forbidden with 68% confidence. If nt = 0, then we measure 
r « 0.3 ± 0.1. If we consider the inflat ionary consistency 
relation, nt — — r/8 (Peiris et al 
constraint r = 0.3l 



-0.1 



(2003)), we then get the 



The uncertainty is smaller, but not 
by much. We predict a 3a detection of PGWs by the zero 
multipole method. 

The detectability of £o using the ideal experiment is 
shown in Fig. ^ For nt = the effective number of a detec- 
tion is cr ~ lOr. We make this approximation by determining 
the detectability for several values of r and then approximat- 
ing a line. For comparison the results are also shown for the 
zero-multipole method with the realistic experiment and for 
measurements of the BB power spectrum with the realis- 
tic experiment described above. We assume we can make 



Figure 7. This is a plot of the distribution of the number of 
positive multipoles for the Monte Carlo simulation for the ideal 
experiment (upper left), the realistic experiment (upper right), 
Planck (lower left), and WMAP (lower right). The dotted red 
line shows where A^-i- = ^Niji„g 



measurements over a range of 60 multipoles for BB mea- 
surements. 

The Monte Carlo simulation for the Wiener filtering 
gives an average of 19 measured TE power spectrum mul- 
tipoles greater than zero out of a total of 52 independent 
multipoles. If the null hypothesis was true, the sign test 
would indicate there is a 3.5% chance of measuring ^ 19 
positive multipoles. This is equivalent to a ~ l.So" detec- 
tion. A plot of the distribution of the number of positive 
multipoles is shown in the upper panel plot of Fig. ^. There 
is an 81% chance for the observed A''+ to give a Icr detection 
of PGWs. 

The S/N test gives a mean value of S/N = —17.1 and 
standard deviation of 7.21. The upper left panel in Fig. ^ 
shows the distribution of the S/N values for the Monte Carlo 
simulation with r = 0.3. If r = 0.3 we would have a 0.8% 
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Figure 8. The S/N statistic distribution for the ideal experiment 
(upper left), realistic experiment (upper right), Planck (lower 
left), and WMAP (lower right). The dotted red line shows where 
S/N = 0. 



Figure 10. This is the plot of the distribution of U for the 
ideal experiment (upper left), realistic experiment (upper right), 
Planck (lower left), and WMAP (lower right). The red dotted line 
is the value for my and the light blue dashed lines enclose the la 
region for U assuming the hypothesis that r = 



I 




Figure 9. This is a plot of (S/N) and crg/ff as a function of r 
for the ideal experiment. The black line is {S/N) and the red line 

is O-g/jv- 



probability of the measured S/N > 0. This negative value 
signifies that a non-zero tensor-to-scalar ratio produced an 
anti-correlation. We can assume that the standard deviation 
would be the same if the mean of S/N was (equivalent to 
r — 0.0), because it is equivalent to adding a constant value 
to every measured value (and hence adding a constant to 
S/N which would not change the error). Therefore, if r = 0, 
the probability of getting S/N < —17.4 is 0.8%, and hence 
we have a 99% chance that r 7^ 0. A plot of (S/N) and o-g/jv 
as a function of r is shown in Fig. ^. As can be seen from 
the plot, we can predict a value of r for any value of S/N. 
The value of (jg/]v is a relatively constant function of r and 
so our prediction about the distribution of S/N for different 
value of r is a good approximation to the true distribution. 

The Wilcoxon rank sum test gives Uavg — rnu = 
— 1.23cri7- The variable Uavg is the mean value for U in the 
Monte Carlo simulations described earlier. The values mu 



o 
Co 
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and ay are given in Section 4.3, The distribution of I] for 



Figure 11. This is the plot of the signal-to-noise ratio (number 
of (Ts) for different values of r for the three different tests. The 
black line is the S/N test, the dashed dark blue line is the sign 
test, and the dotted-dashed light blue line is the Wilcoxon rank 
sum test 



the Monte Carlo simulations with r — 0.3 is shown in Fig. 
p^ . The standard deviation of the distribution of measured 
U is the same as the standard deviation of the distribu- 
tion of [7 assuming the hypothesis that r = 0. The only 
difference between the distributions is that mjj is shifted 
by a constant value. Therefore, there is a 22% chance that 
U — rrijj < —2au- There is also a 40% chance that we mea- 
sure U — mu < —lau, and are not even able to make a la 
detection of PGWs. 

A comparison of the three tests is shown in Fig. 
This is obtained by simulated with with several values of 
r and then interpolating between them. A 2a detection is 
obtained for r = 0.26 {S/N test), r = 0.3 (sign test), and 
r — 0.5 (Wilcoxon rank sum test), highlighting its intended 
use as a monitor of a false positive detection for large r. 
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Figure 12. This figure is tiic same as Fig. ^ but presents for the 
realistic experiment with r = 0.9. 
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Figure 13. An upper limit on r due to the realistic experiment 
with Aio = 10 



5.2 Realistic Ground Based Experiment 

A plot of the error bars for the realistic experiment is shown 
in Fig. |l^ with r = 0.9. Observations on an incomplete sky 
require the multipoles to be binned in sizes of = 10. This 
experiment has much larger error bars than the ideal exper- 
iment and it is not able to detect low values of r with the 
TE cross-correlation only. Plots of the TE power spectrum 
due to density perturbations and PGWs are shown in Fig. 
along with the combined TE power spectrum. 
For this experiment, the constraints on measuring 
are significantly larger than those for the ideal experiment. 
The Icr uncertainty on to is A^o « 10. This corresponds to 
a limit of r < 0.9 with 68% confidence. If we want a 2a 
limit, then the constraint expands to r ^ 1.5. If we assume 
the inflationary consistency relation, then this error on 
would correspond to a Icr upper limit of about r < 0.7. Fig. 
[l^ shows the region of r and nt allowed with 68% confidence 
of lo. 

As mentioned earlier, Fig. ^ shows the signal-to-noise 
ratio of the zero multipole method for the realistic exper- 
iment. The measurements of the BB power spectrum are 



much more sensitive to PGWs and the sensitivity is roughly 
the same as in the ideal experiment. This shows that the 
zero multipole method is less sensitive to PGWs than mea- 
surements of the BB power spectrum. 

The results for the Wiener filtering method are much 
worse than those for the ideal experiment for r = 0.3. Since 
this experiment observes a small portion of the sky, the mul- 
tipoles are correlated and we must bin together to get rea- 
sonably uncorrelated measurements. For this experiment, we 
only have 7 to 8 uncorrelated multipoles, instead of 52 uncor- 
related multipoles in the case where the full sky is observed. 
Getting 7 out of 8 negative multipoles is a 3% probability if 
there are no PGWs. For the Monte Carlo simulations of the 
realistic experiment, on average, half of measured multipoles 
are positive and half are negative. A plot of the distribution 
of the number of positive multipoles is shown in the upper 
right panel of Fig. ^. In this case, we cannot distinguish 
r — 0.3 from r = 0.0 with any significance. 

The S/N test gives an average value of S/N = —0.95 
with standard deviation of 2.64. For the realistic toy exper- 
iment, the distribution of S/N for r — 0.3 is shown in the 
upper right panel of Fig. ^. In order to obtain 68% confidence 
detection of PGWs, we must use r ~ 0.7. In this sense the 
TE test provides monitoring and insurance against false pos- 
itive detection with r > 0.7, which could arise, for example, 
if foregrounds or other systematic effects arer improperly 
removed. 

The last statistical test, the Wilcoxon rank sum test, 
gives Uavg — mu = — 0.20(7 [/. The distribution of U for 
r — 0.3 is shown in the upper right panel of Fig. |l^. This 
gives the weakest result in terms of the three tests for the 
Wiener filtered data. The realistic experiment will not be 
able to constrain r < 0.3 using the TE cross correlation 
power spectrum. Its limit is closer to r < 0.7 — 0.9 at only 
68% confidence depending on the test used. For a higher 
confidence in a detection of PGWs, the value of r would 
need to be much higher. Since the observed distribution of 
U corresponds almost exactly to the simulated distribution 
of U under the assumption that r — 0, therefore we have a 
16% chance of measuring U — mu < —lau- 

5.3 WMAP 

A constraint on r using a measurement of io for WMAP is 
almost impossible. Using error bars consistent with WMAP 
noise, we get A£o « 15 for an input of r = 0.3 and nt — 0. 
The published results of WMAP give limits of r < 0.3 so 
adding this method to the WMAP results would not change 
constraints significantly. In fact, using the real WMAP datan 
we get £o ~ 48. With an uncertainty of Aio « 15, the 
probability of getting a value farther away from £o = 53 
is larger than 50%, so we cannot detect primordial gravita- 
tional waves in the published WMAP data using the zero 
multipole method. 

The results of the Wiener filtering showed that the 
WMAP cannot make a detection of gravitational waves us- 
ing the TE cross correlation power spectrum alone. As with 
the two toy experiments, the result of the scalar and tensor 
separation was similar. The Monte Carlo simulation gave on 

http://lambda.gsfc.nasa.gov/ 
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average gave 13 positive multipoles out of a total of 26 un- 
correlated multipoles. We would get the same result if the 
input data had r — 0.0 so we cannot detect PGWs with 
WMAP using only the TE power spectrum. A plot of the 
distribution of the number of positive multipoles is shown 
in the lower right panel of Fig. ^ As can be seen, this dis- 
tribution of A*'+ for WMAP noise and r = 0.3 is simply the 
distribution for r = 0. 

For WMAP, the S/N test gives the value of S/N = 
—0.02 with a standard deviation of 5.09. The distribution is 
shown in the lower right panel of Fig. ^ The distribution 
is centered around S/N = so there is no chance of using 
this test to detect PGWs in WMAP's TE power spectrum. 
The probability of getting a la or 2a detection is the same 
probability that we would randomly get a detection if there 
are no PGWs. 

The rank sum test gives a value of Uavg — rnu = 
—0.004(7(7, which is implies no ability to distinguish 
WMAP's observed TE data from a data set with no PGWs. 
A plot of the distribution of U for WMAP error bars is shown 
in the lower right panel of Fig. |l^. We reach the same conclu- 
sion for WMAP noise as for the realistic experiment. There 
is only a 16% chance that we can measure U — mu < —lau 
and make a la detection of r = 0.3 

The published WMAP results show an anti-correlation 
of TE power spectrum at large scales. Un fortunately this 



is not a detection of PGWs as theorized in Baskaran et al 



(2006). The contribution to the TE power spectrum due 
to PGWs only changes sign once for £ < 90. If a claimed 
evidence for gravitational waves is to be believed, then the 
TE power spectrum would have to change sign three times 
for £ < 60. In fact, other than the two ant icorrelat ions at 
low £, the rest of the multipoles, up to ^ = 53, are consistent 
with r = 0. None of the described tests applied to the current 
WMAP data will give any detection of PGWs. 



U — mu < —lau and a 9% probability that we measure 
U — mu < —2au for Planck. 



6 COMPARISON OF MEASUREMENTS OF 
THE TE POWER SPECTRUM WITH THE 
BB POWER SPECTRUM 

As mentioned earlier, it was originally suggested that it 
might be easier to detect PGWs using the TE power spec- 
trum instead of the BB power spectrum. For both methods, 
this turned out not to be true. The reason for this is be- 
cause we are trying to measure the TE power spectrum at 
the place where the signal is lowest {Cj^ = 0). In measure- 
ments of the BB power spectrum, if we neglect instrumental 
noise, the signal decreases with a decrease in r and so does 
the cosmic variance limited uncertainty. This is not the case 
for the TE power spectrum. The uncertainty in the mea- 
surement of the TE power spectrum due to PGWs is deter- 
mined by the total TE, TT, and EE power spectra. When 
the TE power spectrum goes to zero, the TT and EE power 
spectrum do not approach zero (in fact, they increase as we 
approach to to). We therefore have a low signal-to-noise ra- 
tio around £o making it very hard to detect PGWs using the 
zero multipole method. Below we give simple summarizing 
arguments why the same is true for the Wiener filtering of 
the TE power spectrum 

If Ni <^ Cf^, the signal-to-noise ratio for the BB power 
spectrum is 



{S/N)i 



= 7 



(20) 



where 



5.4 Planck 

The uncertainty in £o is much better for Planck than for 
the realistic experiment and about twice as large for the 
ideal experiment. The Monte Carlo simulations resulted in 
A^o ~ 3.75 for an input TE power spectrum with r — 0.3 
and nt = 0. This results in « 68% confidence that r 7^ 0, 
under the assumption that nt — 0. 

The sign test gives on average 10 positive measurements 
of the TE power spectrum out of a total of 26 uncorrelated 
multipoles. There is a 16% chance of getting ^ 10 positive 
multipoles if r = 0. A plot of the distribution of the number 
of positive multipoles for Planck is shown in the lower left 
panel of Fig. |^. There is a 50% chance that we will measure 
A''+ < 10 and hence have a la detection of r = 0.3. 

The S/N test gives a value of S/N = -6.24 with a 
standard deviation of 5.09. There is only a 10% chance that 
the S/N test results in a value of S/N larger than zero, if r = 
0.3, and a 10% chance getting S/N < -3.12 if r = 0. This 
is close to a 90% probability of detection. The distribution 
of the S/N variable is shown in lower left panel of Fig. |^. 

Again, the rank sum test gives the lowest confidence 
result with a value of Uavg — mu ~ —0.66au- A plot of 
the distribution of U is shown in the lower left panel of 
Fig. There is a 37% probability that we will measure 
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{2£+l)Uy 



(21) 



If Ne > Cf^ then we wiU not be able to detect PGWs and a 
comparison with the TE power spectrum is not worthwhile. 

If Ni < Cf ^ and r < 1, for the TE power spectrum, 
the signal-to-noise ratio is 



{S/Nh 



\/27- 



[{CJ^Y + + Nt/2){Cf^ + Ni\ 



,1/2 



%/2- 



V2- 



a + f3r 



(22) 



where a and /3 are 



r)TE 

TE riTE 
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where 



^^XY r^XY / 

Di = Ci,t /r 



(24) 



One can see that a and (3 are on the order of unity. There- 
fore, the signal-to-noise ratio is approximated as 



{S/N)te = ^7 



^7- 



(25) 



In other words ifr<a//3~l, BB measurements have the 
obvious advantage in comparison with the Wiener filtering 
of the TE power spectrum. Indeed if r < 1, {S/N)bb ~ 7, 
while {S/N)te ~ 7r < 7. This is because in BB measure- 
ments, applying proper data analysis, we can entirely elim- 
inate contributions of scalar perturbations to CMB polar- 
ization signal as well as to the uncertainties. For the perfect 
Wiener filtering of the TE power spectrum, we can eliminate 
the contribution of scalar perturbations to the signal only, 
but cannot eliminate their contribution to the uncertainties. 



7 CONCLUSION 

The measurement of where the TE cross correlation first 
changes sign can be used to detect or put constraints on 
PGWs. Such constraints are not as strong as the ones given 
by measurements of the BB power spectrum, however it is 
useful to have a supplementary method to detect PGWs. 
We have shown how well the TE mode can constrain the 
amount of PGWs from just a measurement of the angular 
scale where it first changes sign for two different toy experi- 
ments and two real satellite experiments. The absolute best 
limit with which we can measure £0 only gives us less than a 
3(T detection of the PGW component if r = 0.3. The current 
confidence limits gives us r < 0.3 at 95% confidence level. 
Current and future experiments are optimized to measure 
the BB power spectrum if r 0.1 even in the presence of 
foregrounds, which are not taken into account in this pa- 
per. Future satellite experiments should be able to detect 
r < 0.01 which is 10 times better than the sensitivity to r 
than the result of the ideal experiment. If one neglects even 
cosmic variance, the discreteness of £ limits the calculation 
of £0, and the sensitivity to r, to values considerably larger 
than 0.01. The cosmic variance is largest at low £ and is pro- 
portional to the total power spectrum. Since the TE cross 
correlation has contributions from density perturbations the 
errors in the measured TE power spectrum make detecting 
deviations of £0 from 53 difficult, though they also provide 
insurance against a false detection or imperfect subtraction 
of instrumental and foreground systematic effects. 

The other method described in this paper is one in 
which we filter out the signal due to density perturbations, 
leaving only the contribution to the TE power spectrum due 
to PGWs. We then test the resulting TE power spectrum to 
see if it is negative. Three different statistical tests were used 
to see if there was a significant detection of PGWs. The S/N 
test can give a value for r using a comparison with Monte 
Carlo simulations, while the Wilcoxon rank sum test can 
only give an allowable range for r. The sign test will only 
tell us if r 7^ 0. 

Using the Wiener filtering method, we are unable to 
make as significant of a detection as using the zero multipole 



method. The best result was for the S/N test which would 
give a 2.3cr detection of r = 0.3. To detect PGWs on the 
level of 3(7, the tensor-to-scalae ratio r should be r ^ 0.4. 
The sign test would give 2a detection for r = 0.3 and a 
3(T detection for r = 0.45. The Wilcoxon ranked sum test 
gives only a 1.2a detection for r = 0.3 and a Sa detection 
for r = 0.7. Similar results were gotten for the other three 
experiments tested. Thus in the sense of potential to detect 
PGWs, the zero multipole method is the best, next best 
is the S/N test, then the sign test, and the worst is the 
Wilcoxon ranked sum test. 



Baskaran et al. (2006) present illustrative examples in 
which high r is consistent with measured TT, EE, and TE 
correlations. The value of r is so high in these examples that 
if PGWs with such r really existed, current BB experiments 
would already detect PGWs. All models predict that the TE 
cross correlation power spectrum change sign only once for 
£ < 100. The fact WMAP cannot exclude several multipoles 
with Cf^ > in between multipoles of Cf^ < means that 
the TE cross correlation power spectrum either changes sign 
several times for £ < 100 or there is some instrumental noise 
which causes some anticorrelation measurements. Using in- 
strumental noise consistent with WMAP, our Monte Carlo 
simulations give A£o « 16 and £0 > 40, which means that 
there is no evidence of PGWs in the TE correlation power 
spectrum. 
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